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ABSTRACT 



i-^ Context. Longitudinal filament oscillations recently attracted more and more attention, while the 

I ^ restoring force and the damping mechanisms are still elusive. 

< I ' Aims. In this paper, we intend to investigate the underlying physics for coherent longitudinal os- 

-*— > ' 
C/3 ' dilations of the entire filament body, including their triggering mechamsm, dominant restoring 

force, and damping mechanisms. 

Methods. With the MPI-AMRVAC code, we carry out radiative hydrodynamic numerical sim- 

^ . ulations of the longitudinal prominence oscillations. Two types of perturbations, i.e., impulsive 

—^ , heating at one leg of the loop and an impulsive momentum deposition are introduced to the 

f~^ , prominence, which then starts to oscillate. We study the resulting oscillations for a large param- 



eter scan, including the chromospheric heating duration, initial velocity of the prominence, and 
field line geometry. 



o 

CO ' Results. It is found that both microflare-sized impulsive heating at one leg of the loop and a sud- 

denly imposed velocity perturbation can propel the prominence to oscillate along the magnetic 
dip. An extensive parameter survey results in a scaling law, showing that the period of the oscil- 

p^ , lation, which weakly depends on the length and height of the prominence, and the amplitude of 

jrt , the perturbations, scales with -\jR/gQ, where R represents the curvature radius of the dip, and gg 

is the gravitational acceleration of the Sun. This is consistent with the linear theory of a pendu- 
lum, which implies that the field-aligned component of gravity is the main restoring force for the 
prominence longitudinal oscillations, as confirmed by the force analysis. However, the gas pres- 
sure gradient becomes non-negligible for short prominences. The oscillation damps with time in 
the presence of non-adiabatic processes. Compared to heat conduction, the radiative cooling is the 
dominant factor leading to the damping. A scaling law for the damping timescale is derived, i.e., 
T ~ /'■^^£)''^^w"'"^'vg''-^'', showing strong dependence on the prominence length /, the geometry 
of the magnetic dip (characterized by the depth D and the width w), and the velocity perturbation 
amplitude vq. The larger the amplitude, the faster the oscillation damps. It is also found that mass 
drainage significantly reduces the damping timescale when the perturbation is too strong. 
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1. Introduction 

Solar prominences, or filaments when appearing on the solar disk, are cold and dense plasmas sus- 
pended in the corona (Tandberg-Hanssen I 19951 Labrosse et al. 120101 Mackay et al. 120101) . They are 
formed above the magnetic polarity inversion lines. The denser mateiial is believed to be supported 
by the magnetic tension force of the dip-shaped magnetic field lines (Kippenhahn & Schluter|T957l 
Kuperus & Raadu FT974l Guo et al. [20T0t Zhang et al. [20T2t Xu et al. l20T2l Su & van Ballegooijen 
I2012I I. These fascinating phenomena attracted a lot of simulation eff'orts from different aspects, such 
as their formation, oscillations, and eruptions. With respect to the formation, the chromospheric 
evaporation plus coronal condensation model has been studied widely with one-dimensional (ID) 
simulations (e.g., Miiller et al. l2004l Karpen et al. l2005ll2006l Karpen & Antiochos l2008l Antolin 
et al. 120101 Xia et al. 120111 Luna et al. I2012bl l. where no back-reaction on the field topology is 
accounted for It was then for the first time extended to 2.5D by Xia et al. (I2012I I who simulated 
the in situ formation of a filament in a sheared magnetic arcade and showed that the condensation 
self-consistently forms magnetic dips while ensuring force-balance states. This finding strengthens 
the hitherto invariably ID analysis performed for prominence formation and evolutions, as adopted 
by many authors to date. Once a prominence is formed, it might be triggered to deviate from its 
equilibrium position and start to oscillate. 

Observations demonstrate that prominences are hardly static. Besides small-amplitude oscilla- 
tions (Okamoto et al. 120071 Ning et al. 120091 ). large-amplitude and long-period prominence oscil- 
lations have been observed (e.g., Eto et al. l2002l Isobe & Tripathi l2006l Gilbert et al. l2008l Chen 
et al. 120081 Tripathi et al. |2009l Hershaw et al. l20m Bocchiahni et al. lSOTlT) . The observations 
of the prominence oscillations led to the comprehensive research topic of prominence seismology 
(Blokland & Keppens [20nil201 Ibt Arregui & Ballester llom Arregui et al. l20T2l Luna & Karpen 
120121 Luna et al. I2012al l. and the long-term oscillations were considered as one of the precursors 
for coronal mass ejections (CMEs; Chen, Innes, & Solanki i2008i: Chen i2011i) . Of particular interest 
in this paper are the longitudinal oscillations along the axis of prominences/filaments, which were 
first presented in the simulation results of Antiochos et al. (120001 1 discovered from Ha observations 
by Jing et al. (120031) . The phenomenon was further investigated by Jing et al. (I2006I I and Vrsnak et 
al. (120071) . Such large-amplitude oscillations are triggered by small-scale solar eiTiptions near the 
footpoints of the main filaments, such as mini-filament eruptions, subflares, and flares. The initial 
velocities of the oscillations are 30-100 km s ' . The oscillation period ranges from 40 min to 160 
min and the damping times are ~2-5 times the oscillation period (Jing et al. l2006l l. 

Unlike the transverse oscillations whose restoring force is known to be the magnetic tension 
force, the dominant restoring force for the longitudinal oscillations still await to be claiified. Jing et 
al. (I2003I I proposed several candidates for the restoring force, i.e., gravity, the pressure imbalance, 
and the magnetic tension force. Vrsnak et al. (12007) suggested that the restoring force is the mag- 
netic pressure gradient along the filament axis. With radiative hydrodynamic simulations, Luna & 
Karpen (12012!) and Zhang et al. (I2012I I suggested that the gravity component along the magnetic 
field is the main restoring force. Li & Zhang (12012) . on the other hand, suggested that both gravity 
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and magnetic tension force contribute to the restoring force. As for the damping mechanism, it 
really depends on the oscillation mode. For the vertical oscillations, Hyder (1 19661 1 proposed that 
the magnetic viscosity contributes to the decay. For the horizontal transverse oscillations, Kleczek 
& Kuperus d 19691 1 proposed that the induced compressional wave in the surrounding corona acts to 
seemingly dissipate the oscillatory power. More damping mechanisms have been proposed, such 
as thermal conduction, radiation, ion-neutral collisions, resonant absorption, and wave leakage (see 
Arregui et al. 120 121 and Tripathi et al. 12009 1 for reviews). For the longitudinal oscillations, Zhang 
et al. (I2012I I found that non-adiabatic terms such as the radiation and the heat conduction con- 
tribute to the damping, but they might not be sufficient to explain the observed shorter timescale. 
In their simulations the chromospheric heating is switched off, so that the prominence mass was 
nearly fixed. On the contrary, Luna & Karpen (12012) studied the prominence oscillations while 
keeping the chromospheric heating and the resulting chromospheric evaporation. As a result, the 
prominence was growing in length and mass during oscillations. They found that there are two 
damping timescales, a short one for the initial stage and a longer one later The analytical solution 
indicates that the mass accumulation can explain the fast damping of the initial state. As for the 
later slower damping, they suggested non-adiabatic effects such as radiation and heat conduction. 
A quantitative survey is in order to clarify how different geometrical and physical parameters of 
the prominence affect the damping timescale. 

Within the framework of gravity serving as the restoring force for the filament longitudinal 
oscillations, in this paper we try to do a parameter survey, aiming to clarify how the geometry 
of the magnetic field affects the oscillation period and how the combined effects of radiation and 
heat conduction contribute to the damping of the oscillations. We describe the numerical method 
in Section|2] After showing the effects of the perturbation type in Section|3] we display the results 
of our parameter survey in Section^ Discussions and summary are presented in Sections |5]and|6] 

2. Numerical method 

High-resolution observations indicate that a filament/prominence is made of many thin threads 
which are believed to be aligned to the individual magnetic tubes (Lin et al. 120051 1. Since the 
magnetic field inside the filament is quite strong (Schmieder & Aulanier l2012l l and the plasma beta 
is veiy low (y? ~ 0.01 - 0.1) (Antiochos et al. 2000; DeVore & Antiochos l2000l Aulanier et al. 
I2OO6I 1. plus that the thermal conduction is strongly prevented across the field lines, the dynamics 
inside different magnetic tubes can be considered to be independent. Therefore, the formation and 
evolution of a filament thread can be treated as a ID hydrodynamic problem. Following Xia et 
al. (1201 11 1, the ID radiative hydrodynamic equations, shown as follows, are numerically solved by 
the state-of-the-art MPI-Adaptive Mesh Refinement- Versatile Advection Code (MPI-AMRVAC; 
Keppens et al. l2003ll20T2b . 

Jj(pv) + j-^ipv^ + p) ^ pg\\{s) , (2) 

T + T^^^ + ^^) = PS\\v + ^ ~ nnti^MT) + -r-iK-r-) , (3) 

at OS OS OS 
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where p is the mass density, T is the temperature, s is the distance along the loop, v is the velocity 
of plasma, p is the gas pressure, e - pv'^/2 + pfiy - 1) is the total energy density, hh is the number 
density of hydrogen, Wg is the number density of electrons, and g\\{s) is the component of gravity at 
a distance s along the magnetic loop, which is determined by the geometry of the magnetic loop. 
Furthermore, y = 5/3 is the ratio of the specific heats, A(r) is the radiative loss coefficient for the 
optically thin emission, H{s) is the volumetric heating rate, and k - lO^^T^^^ ergs cm ' s ' K ' 
is the Spitzer heat conductivity. As done in previous works mentioned in ^ we assume a fully 
ionized plasma and adopt the one-fluid model. Considering the helium abundance (nHe/ne = 0.1), 
we take p - lAm^nn and p - 23n^k^T, where rrip is the proton mass and k^ is the Boltzmann 
constant. Note that the above equations are difl'erent from those in Luna & Karpen (120121 1 in that 
a uniform cross section is assumed here for the flux tube for simplicity, where expanding flux 
tubes based on given, immobile 3D magnetic fields are adopted in Luna & Karpen ( I2012I I. The 
radiative hydrodynamic equations ([TJQ are numerically solved by the MPI-AMRVAC code, where 
the heat conduction term is solved with an implicit scheme separately from other terms. To include 
the radiative loss, we take the second-order polynomial interpolation to compile a high resolution 
table based on the radiative loss calculations using updated element abundances and better atomic 
models over a wide temperature range (Colgan et al. I2008t . The corresponding values in this table 
are systematically ~2 times larger than the previous radiative loss function adopted by Luna & 
Karpen ([20T2b . 

It is often believed that a prominence is hosted at the dip of a magnetic loop, supported by 
the magnetic tension force. Therefore, we adopt a loop geometry with a magnetic dip, which is 
symmetric about the midpoint, as shown in Fig. [T] The loop consists of two vertical legs with a 
length of s\ , two quarter-circular shoulders with a radius r (the length of each arc, ^2 - ■si , is r:r/2), 
and a quasi-sinusoidal-shaped dip with a half-length of w. The height of the dip is expressed as 
y — D - D cos(7rx/2w) if the local coordinates {x, y) are centered at the midpoint of the dip. The 
dip has a depth of D below the apex of the loop. Such a geometry determines the field-aligned 
component of the gravity, whose distribution along the left half of the magnetic loop is expressed 
as follows: 



^ll(«) = 



In s - Si \ ^ . 

-^0COs(2sjTrs7J, si<s<S2; (4) 



where the gravity at the solar surface g^ - 2.7 x 10^ m s"^, the total length of the loop L, the 
length of each vertical segment ii = 5 Mm, and S2-s\+ nrll Mm. The total length of the dip is 
2w - L — 2s2- The field-aligned component of the gravity in the right half is symmetric to the left 
half. The parameter h - s\ + r — D gives the height of the central dip above the lower boundary. 

Our simulations start from a thermal and force-balanced equiUbrium state where the back- 
ground heating is balanced by radiative loss and thermal conduction, and the plasma in the loop is 
quiescent. The simulations are divided into three steps: (1) Prominence formation: A prominence 
forms and grows near the center of the magnetic dip as chromospheric material is evaporated into 
the corona and condensates due to thermal instability after chromospheric heating is introduced 
near the footpoints of the loop; (2) Prominence relaxation: The prominence relaxes to a thermal 
and force-balanced equilibrium state as the localized heating is halted and the chromospheric evap- 
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Fig. 1. Magnetic loop used for the ID radiative hydrodynamic simulations of the prominence os- 
cillations. Note that the horizontal and the vertical sizes are not to scale. 



oration ceases; (3) Prominence oscillation subjected to perturbations: The prominence starts to 
oscillate with a damping amplitude after perturbations are introduced. In step 1, which lasts for a 
time interval of Afi, the heating term H(s) in Eq. Q is composed of two terms, i.e., the steady 
background heating Hq{s) and the localized chromospheric heating Hi{s), which are expressed as 
foUows: 



Hois) 



Eo exp{-s/H„), 



s < L/2; 



Eo exp[-(L - s)/H^l L/2 ^ s < L; 



(5) 



Hi(s) 



El, s < s,r; 

El exp[-(i - s,r)/A], s,r < s ^ L/2; 

El exp[-(L - s,r - s)/A], L/2 < s ^ L- s,. 

El, S > L - Str', 



(6) 



where the quiescent heating term Ho is adopted to maintain the hot corona with the amplitude 
£■() = 3 X 10 '^ ergs cm"^ s ' and the scale-height H„, - L/2, and the locaUzed heating term Hi is 
adopted to generate chromospheric evaporation into the corona with the amplitude Ei = 10"^ ergs 
cm"-' s"', the transition region height s,r = 6 Mm, and the scale height A - 10 Mm. The heating 
is taken to be symmetric in order to form a static prominence near the magnetic dip center, so that 
we can easily control the manner how the prominence is triggered to oscillate. Our methodology is 
different from Luna & Karpen (I2012I I who used asymmetric heating which spontaneously leads to 
the oscillation once the prominence is formed. In step 2, Hi is switched off. Owing to the absence 
of the chromospheric evaporation, the gas pressure inside the magnetic loop drops down, so the 
compressed prominence expands until a new equilibrium is reached, which roughly takes less than 
2.4 hr. In step 3, a perturbation is introduced to the prominence in order to trigger its oscillation. 
Note that Ho remains throughout the simulations. 

From the observational point of view, there might be two kinds of perturbations. The first one 
is an impulsive momentum injected to the magnetic loop as the magnetic reconnection near the 
footpoints rearranges the magnetic loop rapidly. The second is impulsive heating due to subflares 
(e.g., Jing et al. l2003l Vrsnak et al. [20071 Li & Zhang |2012] | or microflares (Fang et al. 12006^ 
near the footpoints of the magnetic loop where a large amount of magnetic energy is impulsively 
released through magnetic reconnection. The gas pressure is greatly increased that could propel the 
prominence to oscillate along the dip-shaped field lines. In our ID simulations, we separate the two 
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effects to see their difference. In one case, a velocity perturbation with the following distribution is 
imposed to the prominence, 

0, i < ipi - 6; 

vo(s - Spi + 6)16, Sp\- 5 < s ^ ipi; 
v{s) = \ vo, Spi < 5 < v; (7) 

Vq{-S + Spr + 5)15, ipr < 5 < ipr + 6; 

0, s > ip,-, 

where Spi and s^^ are the coordinates of the left and right boundaries of the prominence, 5=10 
is the buffer zone which allows that the perturbation velocity varies smoothly in space, and vq is 
the perturbation amplitude. In the other case, impulsive heating {H2), as described as follows, is 
introduced near the right-hand footpoint of the magnetic loop, 

W ^peak) \} ^peak) 



H2{s) = £2 exp 



scale scale 



(8) 



where the heating spatial scale Sscaie - 2.5 Mm, the peak location Sp^ak - 245 Mm, the heating 
timescale t scale = 5 min, and the peak time tpgak = 15 min. The heating ramps up to the peak for 15 
min and then fades down to 0. 

As for the boundary conditions, all variables at the two footpoints of the magnetic loop are fixed, 
which is justified because the density in the low atmosphere is more than four orders of magnitude 
higher than that in the corona. The same approach has been adopted by Ofman & Wang (I2002I I and 
Xia et al. (1201 11 1, assuming that the coronal dynamics has little effect on the low atmosphere. The 
approach was verified by Hood (11986) with the parameters being far from the marginal stability. 
The violation of the rigid wall conditions in certain cases was discussed by van der Linden et al. 
(fT994l i. 

3. Effects of the perturbation type 

In order to check how the two types of perturbations as described in ^influence the characteristics 
of the prominence oscillations, we perform simulations of the oscillations which are excited by the 
two types of perturbations while keeping Afi - 7.2 hr, r = 20 Mm, D = 10 Mm, and L = 260 Mm. 
In case A, the prominence oscillation is triggered by a velocity perturbation over the whole 
prominence body. With vq - -40 km s ' (the minus means that the velocity is toward the left), the 
temporal evolution of the plasma temperature distribution along the magnetic loop is displayed in 
the left panel of Fig. |2] It is seen that in response to the perturbation, the prominence, signified by 
the low temperature, starts to oscillate around the equilibrium position. The oscillation amplitude 
decays with time. Fitting the trajectory of the mass center of the oscillating prominence with a 
decayed sine function 

2n 
s ^ so+Aosm(—t + 4>o)exp(-t/T), (9) 

we find the initial amplitude Aq - 34.9 Mm, the oscillation period P - 84.3 min, and the damp- 
ing timescale r - 272 min. Assuming that the prominence thread has a cross-section area of 
~ 3.14 X 10'"^ cm^ (Lin et al. i2005i) . the initial kinetic energy of the oscillating prominence thread 
is estimated to be ~ 7.2 x 10^^ ergs. It is noted that the single decayed sine function, as used for 
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Fig. 2. Comparison of the evolutions of the temperature of the loop between the two types of 
perturbations. The left panel corresponds to the case with velocity perturbations with vo = -40 
km s ' and the right panel to the case with localized heating perturbations with £2 = 0.24 ergs 



cm ^ s ' . 



fitting the Ha observations (Jing et al. 120031 Vrsnak et al. 120071 Zhang et al. 120121 ). fits the simu- 
lated observations very well. On the contrary, a combination of Bessel function and an exponential 
decay function is necessary to fit the initial overtone in the simulations of Luna & Karpen (120121) . 
which results from the continual mass accumulation. 

In case B, the prominence oscillation is triggered by the impulsive heating which is deposited 
near the right leg of the magnetic loop in order to mimic a microflare near the prominence. To do 
that, an impulsive heating term H2(s) in Eq. ^ is added to the heating term H in Eq. (O, where 
Speak - 245 Mm meaning the heating is concentrated at a height of 15 Mm above the right footpoint 
of the magnetic loop. 

The right panel of Fig. |2]depicts the temporal evolution of the temperature distribution along the 
magnetic loop with E2 = 0.24 ergs cm"^ s '. With the typical cross-section area of a prominence 
thread being ~3. 14x10'"* cm^, the corresponding total energy deposited into the single magnetic 
loop Eheating IS 1.8x10^^ crgs. This value is reasonable since observations indicate that the total 
energy of a microflare is 10^^-10^^ ergs or even more (e.g., Shimizu et al. 120021 Hannah et al. 
120081 Fang et al. 1 20 101 ). and several percent of the released energy goes into one prominence thread. 
From another point of view, under the framework of magnetic reconnection model for microflares, 
the magnetic energy release rate is estimated to be B^Vi„l{AnL). With the magnetic field B ~ 20 
G, the reconnection inflow speed v,„ being about 0.1 times the Alfven speed which is about 1000 
km s ' (Jiang et al. 120121 ). and the spatial size L - 10", the energy release rate is estimated to be 
~ 0.88 ergs cm"^ s ', which is of the order adopted here. Fitting the trajectory of the oscillating 
prominence with the damped sine function as shown in Eq. (|9]) yields A^) - 35.8 Mm, P - 84.3 
min, and r - 268 min. The corresponding initial velocity is also -40 km s '. This indicates that 
a typical microflare near the leg of the magnetic loop hosting a prominence thread can excite the 
prominence longitudinal oscillations with an initial velocity of tens of km s '. The corresponding 
kinetic energy is only ~ 7.2 x 10^^/1.8 x 10^^, i.e., ~4% of the deposited thermal energy. The 
remaining ~96% of the energy deposit contributes to the heating of the chromosphere. 
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4. Parameter survey 

The results in ^reveal that the oscillation period does not strongly depend on the two types of 
perturbations, i.e., impulsive momentum and localized heating at one footpoint used in our inves- 
tigation. Note that we concentrate on the oscillation characteristics which follow the small tran- 
sient/excitation phase, already obtained from simple decaying sinusoidal fitting. A small difference 
in the decay timescale exists between the two perturbation types. With the same initial velocity, 
the decay timescale is 4 minutes shorter in the case of impulsive heating than that in the case of 
impulsive momentum. However, the relative variation, 1.4%, is very small. Therefore, we can con- 
clude that the oscillation is basically intrinsic and the characteristics of the oscillation depend on 
the prominence itself and the geometry of the magnetic loop in our case where there is no mass 
accumulation and the oscillations are excited by either impulsive momentum or localized heating. 
The prominence feature is only characterized by the thread length (0, and the geometry of the 
magnetic loop is characterized by r, D, and w as depicted in Fig. [T| Among the three geometrical 
parameters, h - s\ + r- D determines the height of the prominence, D and w determine the curva- 
ture of the magnetic dip. If other parameters are fixed, the length of the prominence is determined 
by the duration of the chromospheric evaporation in step 1, i.e., Afi, as described in ^ Besides, 
the decay timescale might vary with the perturbation amplitude, therefore another parameter is the 
initial perturbation velocity vq. In this section, we perform a parameter survey to investigate how 
each individual one among the five parameters (Af i , r, D, w, and vq) changes the oscillation period 
and the decay timescale. For each parameter, several cases with different values are simulated with 
other parameters fixed. In our simulations, we set r = 10 Mm, D - 5 Mm, w = 110 Mm, and 
vo = -20 km s"' when varying Afi. We set Afi = 7.16 hr, D = 5 Mm, w = 90 Mm, and vq - -20 
km s"' when varying r. We set Afi - 7.16 hr, D = 5 Mm, r = 10 Mm, and vq = -20 km s"' when 
varying w. We set Afi = 7.16 hr, r = 20 Mm, w - 93.6 Mm, and vq - -20 km s"' when varying 
D. We set Afi = 7.16 hr, r = 20 Mm, w = 93.6 Mm, and D = 10 Mm when varying vq. Since the 
oscillation characteristics are found nearly insensitive to the perturbation type, we use the velocity 
perturbation to excite the oscillations in the survey. 

4.1. Length and mass of the prominence 

After finishing the first two steps of the simulations as described in ^ we get a quasi-static promi- 
nence. The dependence of the prominence length I on Afi , h, D, and w is shown in the four panels 
of the upper row of Fig. |3] It can be seen that I, which fits into the scaling law / ~ Af j'^", increases 
with the duration of the heating time Af i . It is understandable since more chromospheric plasma is 
evaporated into the corona when Afi increases. The length I decreases with has I ~ /z "^^, which is 
probably because it takes a longer time for the more tenuous corona to condensate as the height of 
the magnetic dip increases, and therefore the effective heating time is shorter The length / decreases 
with D as I ~ D "^', which can be understood as the prominence becomes more compressed as the 
magnetic dip becomes deeper However, the length of the prominence does not vary considerably 
with w. Of course, w should not be too small, otherwise thermal instability would not occur The 
lengths of these simulated prominence threads are consistent with the reported values, i.e., tens of 
Mm(Linetal.|2005ll. 
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Fig. 3. Scatter-plots of the total length I (upper panels) and mass M (lower panels) of the promi- 
nences at the end of relaxation step as functions of Af i , h, D, and w. 

The dependence of the prominence mass M on Af i , h, D, and w is shown in the four panels of 
the lower row of Fig. |3] It can be seen that the dependence of M on Afi, h, and w is similar to I. 
Their difference is that I decreases with D whereas M does not change with D, which means that 
the plasma number density (10"^- 10" cm"^, and the corresponding density is 10"''*-10"'^ g cm"^) 
is higher in the prominence with a deeper magnetic dip. A scaling law is obtained by fitting the 
data points, which is M ~ Af°^^/i"°-''^. 

It is noted that the above results are derived for a dipped magnetic loop filled via chromospheric 
evaporation with a limited lifetime, where the prominence thread can sustain in the corona. In the 
case of magnetic loops without a dip (e.g., Mendoza-Briceno et al. l2005l l or with a shallow dip and 
asymmetric heating (e.g., Karpen et al. I2006I I. condensations repetitively form, stream along the 
magnetic field, and ultimately disappear after falling back to the nearest footpoint. Therefore, the 
mass and length of the prominence evolve dynamically, without reaching an equilibrium value. 



4.2. The oscillation period and decay iimescale 

As the velocity perturbation is introduced to the quasi-static prominence, the prominence starts 
to oscillate. Fitting the trajectory of the oscillating prominence with the damped sine function 
shown in Eq. (|9|l, we get the oscillation period (P) and the decay timescale (t) for each case in the 
parameter survey. 

The variations of P along with the parameters /, h, D, w, and vq are shown in the upper row of 
Fig. |4] It is seen that P increases slightly with / and vq, and decreases slightly with h. However, it 
increases seriously with w and decreases with D. To fit the variations with a scaling law, we obtain 



f 



i().16i,-0.05r(-0.54,„0.9I,,0.05 



=D 



w-^'v" 



Therefore, the period of prominence longitudinal oscillations relies 
dominantly on the geometry of the dip, especially its curvature. It is noted that the range of P is in 
agreement with the reported values in previous studies (e.g., Jing et al. l2006l l. 

The variations of r along with the five parameters are shown in the lower row of Fig . |4] It is seen 
that T increases significantly with I and D, and decreases with w and vq. It is noted that in the cases 
of |vo| = 70 and 80 km s"\ part of the prominence mass drains down to the chromosphere, which 
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Fig. 4. Scatter-plots of the period P (upper panels) and damping time t (lower panels) of the promi- 
nences in the oscillation step as functions of I, h, D, w, and vq. The values of P and t in the cases 
|vo| = 70 and 80 km s ' that cause mass drainage at the footpoint of the coronal loop are marked 
with triangles in the right panels. 

is why the triangles in the lower-right panel of Fig. |4]do not follow the trend of the data points 
denoted by the diamonds where |v()| < 70 km s"'. The decay timescale does not vary significantly 
with h. To fit the variations with a scaling law, we obtain r ~ ;i-63/^-o.i8£)0.66|^-i. 21^^-0.30^ where the 
cases with prominence drainage are not included in the fitting. The values of t are also in the same 
order of magnitude as the observed ones. 

5. Discussions 

5. 1 . Restoring force 

For an oscillating phenomenon, the most important thing is the determination of the restoring force, 
which directly decides the oscillation period. In our ID hydrodynamic simulations, the only forces 
exerted on the prominence are the gravity and the gas pressure gradient, both are restoring forces 
for the longitudinal oscillations. In order to compare their importance, we calculate the two forces 
in the case with Afi = 7.16 hr, vo = -40 km s'\ r - 20 Mm, D = 10 Mm, and w - 93.6 Mm. 
The two forces are calculated when the prominence is the furthest from the equilibrium position. 
Despite that the plasma in prominences is hundreds of times denser than the ambient corona, it 
is not an ideal rigid body. For oscillations with higher modes as studied by Luna et al. (I2012al l. 
the pressure gradient changes rapidly along the prominence thread. For the fundamental-mode 
oscillations in this paper, the prominence oscillates as a whole and the pressure gradient changes 
slightly along the thread. Therefore, for simplicity, we compare the overall magnitude of the two 
forces by a simple calculation instead of as point-to-point one in the simulations. The integral of 
the gravity force is quantified between the two ends of the prominence, i.e., Fg - [ " p\g\\\ds - 



Jlefl 

Ii f PSb^I sin( '^'-^~/ )\ds, where a unit area is assumed for the cross section. The integral of 

■right . 



pressure gradient force over the prominence is expressed as Fp - J " \dp/ds\ds - \pright - 
pieftl- The left and right boundaries of the prominence are defined to be where the density drops to 
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Fig. 5. Temporal variation of Fg/Fj, when the displacement of the prominence reaches maximum 
during each half-cycle in the case of r = 20 Mm and D = 10 Mm. The velocity perturbation is -40 

km s"'. 

7x10"'"^ g cm"^. Figure |5]displays the temporal evolution of the ratio Fg/Fp, from which it is seen 
that the gravitational force is generally ~10 times larger than the gas pressure gradient force. 

Since the gravity is the dominant restoring force, the overall motion of the prominence can also 
be described for simplicity as 



d X 
M-^ = M^ll = -MgQ^:^ sin(— ), 



TiD 

2w 



,nx, 
w 



(10) 



where x - s - L/2 is the displacement of the prominence from the equilibrium position. It is not 
easy to solve this equation analytically. However, if the oscillation amplitude is much smaller than 
the half width of the whole magnetic dip (w), we get the approximation sin(:7rx/vv) » nx/w. So, the 
above equation is simplified to be 



d X 



nD nx 



^-nr = -MgQ- 



dt 



2w w 



(11) 



with the solution x - Aq sin(^r + (p). The corresponding period is 



P = 




(12) 




Such a period can also be readily obtained if the prominence is taken in analogy to a pendulum 
whose period is 



(13) 



where R is the curvature radius of the dipped magnetic loop. With the shape of the loop being 
y = Z)-Dcos(7rji:/2w), the curvature radius at the loop center is approximated to be/? = 2w^/{D7t^). 
Substituting R into Eq. ( fTlT l, we get P - -\j8w^/(gQD), the same as Eq. (fT2l i. Figure|6]compares the 
oscillation periods obtained from the hydrodynamic simulations (diamonds) and those estimated 
from Eq. (fT2l i [solid line) when the two parameters, D and w, are changed. It is revealed that 
Eq. ( fT2l l is a very good approximation for estimating the period of the prominence longitudinal 
oscillation. Of course, it should be kept in mind that the derivation of Eq. (fT2l i is based on the 
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assumption that the dipped magnetic loop has a sinusoidal shape. More generally, the oscillation 



period is related to the local curvature radius R by the formula P -2n -\/R/gQ, as also demonstrated 
by Luna & Karpen (l20T2b . 

Recently, Luna et al. (I2012al l extended the theoretical analysis of longitudinal prominence os- 
cillations by including the effect of the pressure gradient force. They found that the ultimate fun- 
damental frequency of the oscillations is found from cij^^^^ = llt^ + w^, where cuig and w, stand for 
the gravity-driven and pressure-driven frequencies, respectively. The ratio of the two frequencies 
a>l/a>^ = Riim/R, where R^^ denotes the critical value of the curvature radius (R) of the magnetic 
dip. If 7? «; Rii„, then gravity dominates over pressure in the restoring force of longitudinal oscil- 
lations. They pointed out that the reported values of the curvature are small compared with /?iim, so 
that it is reasonable to ignore the effect of the pressure term in most cases. In our parameter survey, 
Riim - 0.175(L - 1)1 ranges from 760 to 2100 Mm and the ratio R/R\im ranges from 0.1 to 0.5. 
Hence, their theoretical results of gravity being the main restoring force for the fundamental mode 
in this parameter range are thus confirmed by our simulations . 

For a prominence above the solar limb, all the parameters in Eq. (fT2l) can be roughly measured. 
Combined with the results in this paper, the comparison between simulations and observations 
in Zhang et al. ( I2012I I implies that Eq. (fT2l i is a good approximation to estimate the oscillation 
period. For the prominence longitudinal oscillations on the solar disk, i.e., filament longitudinal 
oscillations, only the oscillation period can be unambiguously measured. Eq. (fT2l l then provides 
a diagnostic tool for inferring the geometry of the dipped magnetic loop. Especially, when w can 
be roughly estimated from force-free magnetic extrapolations, the depth of the dip, D, can be 
determined. At least, we can estimate the curvature radius of the dipped magnetic field, R, through 
Eq. ( fTSl ). After the determination of R, Luna & Karpen (120121) further proposed an approximate 
method to estimate the magnetic field in the prominence. 

Besides the dominant dependence on the geometric parameters, the oscillation period also 
weakly changes with the length and the height of the prominence, as well as with the initial velocity. 
These can be understood as follows: (1) Dependence on the prominence length: As the prominence 
thread is shorter, the ratio of the gas pressure gradient to the gravity would increase as indicated by 
our simulations, therefore, the gas pressure gradient would contribute to the restoring force, result- 
ing in a shorter oscillation period; (2) Dependence on the prominence height: As seen from Fig. [3] 
with other parameters the same, a high prominence has a shorter length. Therefore, with the same 
reason as in (1), the oscillation period would be smaller; (3) Dependence on the initial velocity: 
Since sin(;rx/w) is always smaller than nx/w in Eq. (fTOl l. the nonlinear term would naturally lead 
to a long period as the oscillation amplitude increases. 

5.2. Damping mechanisms 

If the energy dissipation terms such as the radiative cooling and the heat conduction are removed 
from Eq. ^, as we did in a test simulation, we found that the prominence oscillation does not 
damp at all. When the two non-adiabatic terms are kept, the prominence oscillation always damps. 
In order to see the importance of the two terms, we calculate the time integrations of radiative loss 
(Es) and thermal conduction (Eq) of the whole system after subtracting the corresponding values 
when the prominence is static at the center of the dip. Here Er and Eq are the integrals of the 
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Fig. 6. Comparison of the periods of the prominence oscillations from simulations (diamonds) and 
theoretical analysis (solid line) as a function of the depth of the magnetic dip D (left panel) and the 
width of the dip w (right panel). Note that both axes are in logarithmic scale. 
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Fig. 7. Temporal variations of Ej^/Ec in the oscillation step in the cases of vq = -40, -50, and -60 
km s"^ 

radiative and the conductive terms in the energy equation Eq. (O, where the integrals are taken 
in the whole corona above the two footpoints. The evolutions of the ratio (EJ^/Ec) in the cases of 
vo - -40, -50, and -60 km s"' are displayed in Fig.|7] It is seen that the ratio is always larger than 
unity. Especially in the early stage of the oscillation when the amplitude is still large, Er is even 
one order of magnitude larger than Eq- It is also revealed that as the initial velocity increases, /sr 
becomes more and more important in most of the lifetime of the oscillation. Our results support 
the conclusions of Terradas et al. (120011 I2005I I that radiative loss is responsible for the damping of 
the slow mode of prominence oscillations in the dip-shaped magnetic configurations, which seems 
to be different from the case of slow-mode waves propagating in the coronal loops where heat 
conduction contributes more to the damping (De Moortel et al. I2002al I2002bl l. 

The role of the radiative cooling can be understood in a simple model as follows: Since there 
are two segments of the corona in the magnetic loop, as the prominence oscillates, one part would 
be attenuated and the other be compressed. Suppose that the total length of the coronal part of the 
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magnetic loop is unity, which includes the part x, which is to the left of the prominence, and the 
other part I - x, which is to the right of the prominence. Hence, the densities of the corona on the 
two sides are proportional to 1 /x and 1/(1 - x), respectively. The total opticaUy-thin radiative loss 
of the coronal part is proportional to x^^ + (1 - xy^, which is the minimum when x = 0.5, i.e., 
when the prominence is situated at the equilibrium position. Whenever the prominence deviates 
from the loop center, the cooling becomes larger, dissipating the kinetic energy of the oscillating 
prominence. The model is best illustrated by the relationship between the damping timescale (r) 
and the initial amplitude of the oscillation, i.e., Aq in Eq. (|9]l. As Aq increases, one of the two 
coronal parts is more severely compressed, so the radiative cooling x^^ + (1 - xy^ deviates further 
away from the minimum value, i.e., it becomes larger As a result, the oscillation decays more 
rapidly. 

Based on the sinusoidal function, Aq cc vqP. Substituting Eq. (fT2l) into it, we get Aq oc vqwD '^^. 
With this, it is not difficult to understand the positive correlation between the decay timescale r and 
D and the negative correlation between t and w as revealed by the lower row of Fig. |4] Along this 
line of thought, the dependence of the decay timescale on the prominence length can be explained 
as follows: As the prominence thread is longer, the coronal part of the magnetic loop, which radiates 
out the thermal energy, is shorter More importantly, the longer thread, with the same initial velocity, 
has a larger kinetic energy. Therefore, it takes a longer time for the compressed coronal part to 
radiate it out. 

It is seen from the first six cases (i.e., |vo| from 10 km s"' to 60 km s"') in the lower-right panel 
of Fig. |4] that the decay timescale decreases with the initial perturbation velocity nearly linearly. 
However, when vq is larger than 70 km s ', part of the prominence would overpass the magnetic 
loop apex and drain down. The critical velocity for the prominence to reach the loop apex can 



be roughly estimated as v„,f, ~ ^J2gQD - 23 yjD/Mm km s"'. Therefore, the value of Vcriti is 
73 km s"' in the case of D =10 Mm. As revealed from our simulations, even when vo=-70 km 
s ', mass drainage already happens, although the amount of the drainage is much less than that in 
the case of vq =-80 km s"'. The temperature evolution along the loop in the case of vq =-80 km 
s"' is presented in Fig. |8] It is seen that part of the prominence falls down to the left leg of loop, 
leading to the drainage of the prominence mass and kinetic energy as well, while the remaining part 
continues to oscillate along the dip. The oscillation period and the decay timescale in the cases with 
mass drainage are marked as triangles in Fig. E] Their periods, ~90.6 min, are slightly below the 
trend defined by other cases without mass drainage (diamonds), which is consistent with the weak 
positive correlation between P and the prominence length /. However, the damping timescales are 
greatly reduced, compared to the trend defined by other cases without mass drainage as seen from 
the lower-right panel of Fig. |4] Such a result, namely that mass drainage would greatly reduce the 
decay timescale, might explain the mismatch between the simulation and the observation of the 
decay of a prominence oscillation reported in Zhang et al. (I20I2I) . 

6. Summary 

In this paper, we carry out ID hydrodynamic simulations of longitudinal prominence oscillations 
using the MPI-AMRVAC code, extending earlier numerical simulations of prominence formation 
(Xia et al. 1201 11 1 and of prominence oscillations (Luna & Karpen l2012t Zhang et al. I2012I I. The 
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Fig. 8. Temporal evolution of the temperature along the magnetic loop when the initial velocity 
perturbation is as large as vq =-80 km s"'. Note that the prominence overpasses the magnetic loop 
apex and drains down to the chromosphere at the left footpoint around f = 0.8 hr. 

simulations are divided into three steps: First, a prominence forms and grows near the center of 
the dip-shaped coronal loop due to chromospheric heating and the subsequent thermal instability. 
Then, it relaxes to a quiescent state after the chromospheric heating is switched off. Subjected 
to two kinds of perturbations that mimic subflares, the prominence starts to oscillate along the 
dip. Within the framework of the evaporation-condensation model, we obtained scaling-laws for 
the prominence length (/) and mass (M), which are expressed as / ~ Afj-^°/z"°-^^D"°-^' and M ~ 
Af^^^/i " ■^'*, where Af i is the time duration of the chromospheric heating and evaporation, h is the 
prominence height, D is the depth of the magnetic dip. It is found that / is insensitive to the half 
length of the magnetic dip (w) once w is large enough, say, 60 Mm; M is insensitive to D and 
w. Both transient heating at one leg of the loop and an impulsive velocity perturbation applied 
to the prominence as a whole are capable of driving a coherent oscillation along the dip. The 
oscillation properties are found insensitive to the perturbation type in the regimes studied. In the 
case of the transient heating, ~4% of the deposited energy is converted into the kinetic energy of 
the prominence. The longitudinal oscillations are sustained mainly by the tangential component of 
gravity, except when the prominence is short and the gas pressure gradient becomes also important. 
Both simulations and linear analysis reveal that the period of oscillation (P) is In yjRjgQ, where R 
denotes the curvature radius of the dip, as also found by Luna & Karpen (I2012I I. Other parameters, 
such as the length and the height of the prominence, as well as the perturbation velocity, also affect 
P, though slightly. The longitudinal oscillations damp in the presence of non-adiabatic effects, i.e., 
radiative loss and thermal conduction (Soler et al. |2009) l, among which the radiative loss plays a 
leading role. With the parameter survey, we obtained a scaling-law for the decay timescale t, which 
is expressed as r ~ /i 63 £)0.66j^-i. 21^^,^-0.30^ where vo is the initial velocity perturbation. We also found 
that prominence mass drainage, once it happens, significantly reduces the decay timescale, which 
may explain the mismatching between the simulations and the observations disclosed by Zhang et 
al. (I2OT2I 1. 

It is worth mentioning the limitation of the applications of the above results. According to this 
paper, the mass of a prominence thread is insensitive to the depth D and the width w of the magnetic 
dip. This is based on the prominence formation directly via chromospheric evaporation with a fixed 
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lifetime Af i . According to Xia et al. (1201 11 1, the prominence would grow via siphon flow even when 
the localized heating is switched off', though the growth speed is much slower Recently, Luna et 
al. (I2012al l pointed out that the restoring force of the longitudinal oscillations depends on the depth 
of the magnetic dip. For shallow dips, gas pressure plays an important role, while gravity is the 
main factor for deep dips. Besides, Li & Zhang (I2012I I suggested that magnetic tension may also 
contribute to the restoring force. As for the damping mechanisms, several other effects might be 
taken into account in the future simulations, such as the wave leakage and plasma viscosity (Ofman 
& Wang |2002l l. However, some will only be quantifiable in true multidimensional configurations, 
e.g. starting from the prominences formed in Xia et al. (1201 11 1. 
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